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Abstract. In this paper we discuss whether thermodynamical concepts and 
in particular the notion of temperature could be relevant for the dynamics of 
granular systems. We briefly review how a temperature-like quantity can be 
defined and measured in granular media in very different regimes, namely the 
glassy-like, the liquid-like and the granular gas. The common denominator will 
be given by the Fluctuation-Dissipation Theorem, whose validity is explored by 
means of both numerical and experimental techniques. It turns out that, although 
a definition of a temperature is possible in all cases, its interpretation is far from 
being obvious. We discuss the possible perspectives both from the theoretical 
and, more importantly, from the experimental point of view. 
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1. Introduction 

Granular matter EI constitutes one of the most famous example of non- 
equilibrium, athermal systems in which the usual temperature does not play any 
role. However, the fact that they involve many particles brings naturally a strong 
motivation to treat them with thermodynamic or thermodynamic-like methods. It is 
therefore of primary importance to determine whether such attempts are feasible, and 
if usual tools of thermodynamics can be generalized. In particular, can concepts like 
entropy or temperature have any meaning or use ? 

A thermodynamic approach is in general justified when one is able to identify a 
distribution that is left invariant by the dynamics {e.g. the microcanonical ensemble), 
and to assume that this distribution will be reached by the system, under suitable 
conditions of 'ergodicity'. Unfortunately, because energy is lost through collisions or 
internal friction, and gained by a non-thermal source (vibrations, tapping, shearing,...), 
the dynamical equations do not leave the microcanonical or any other known ensemble 
invariant. This raises several questions, among which the most important is whether 
it is possible in principle to construct a coherent thermodynamics for these "non- 
thermal" systems. 

The concept of temperature has been in fact widely used in the context of the 
different states of granular media. The most straightforward definition comes from 
the case of dilute, strongly vibrated granular systems, which reach a non-equilibrium 
stationary state: by analogy with molecular g "granular temperature" Tg can 

be defined in terms of the average local kinetic energy per particle, and treated as 
one of the hydrodynamic fields. This approach can a priori be extended to denser, 
liquid-like, strongly vibrated systems. 

In dense, slowly evolving dense granular media, the situation is more complex 
since reaching a stationary state is experimentally or numerically very difficult. These 
systems actually exhibit aging and memory 017]. Analogies with other aging 

systems have also led to the definition of dynamic temperatures as quantifying the 
violation of the equilibrium Fluctuation-Dissipation Theorem (FDT) ^ ^ . 

In this paper, we will briefly review the different cases of dense, liquid-like and gas- 
like granular media, focusing on the notion of temperature as defined in the framework 
of the Fluctuation-Dissipation Theorem. 

The paper is organized along the following lines. In section [21 we review the 
Fluctuation-Dissipation Theorem and its generalizations in order to provide the tools 
that will be used to test its validity in the different regimes. In Section 01 we discuss 
the case of dense granular media. We first present a simple pedagogical 1-d model for 
a granular medium where the notions of entropy and temperature emerge naturally 
and that could be used as a reference, although unrealistic, example. We then review 
the basics of Edwards' approach as well as some results obtained recently to test 
its validity. Next, in section 01 we describe experiments devoted to the test of the 
Fluctuation-Dissipation Theorem in a vibro-fluidized granular medium, probed by 
means of a torsion pendulum. Section [SI is devoted to the discussion of the case of 
granular gases. Also in this case we shall present some procedures to test the validity 
of the Fluctuation-Dissipation Theorem and to define a notion of temperature. Finally 
Section [HI is devoted to the conclusions as well as to drawing some perspectives. 
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2. Fluctuation-Dissipation Theorem and generalizations 

Let us consider an equilibrium system in contact witli a thermostat at temperature 
T. For the purposes of this paper, the Fluctuation-Dissipation Theorem |H] can be 
seen as relating the random diffusion and the mobility of a tracer particle in a gas or 

liquid: one possible version is the Einstein relation {[r{t) — = 2dT ^^^^^\^^^ ''^ , 

where r is the position of the particle and / is a constant perturbing field, and the 
brackets denote average over realizations. 

Since the FDT is a feature of equilibrium systems, nothing guarantees its validity 
in out of equilibrium systems. The situation is even worse in athermal systems, since 
it is not clear which quantity should play the role of the temperature. 

Recent developments in glass theory, especially those related to their out of 
equilibrium dynamics, have shown that in slowly evolving, aging systems, the FDT is 
in fact modified in a very interesting fashion (for a recent review see In a class of 

mean-field models, which contains, although in a rather schematic way, the essentials 
of glassy phenomena |11[I12| . and whose aging dynamics was solved analytically |13j . 
a feature that emerged was the existence of a temperature T^yn for all the slow modes 
(corresponding to structural rearrangements) |14l I15| . This dynamical temperature 
T^yn in fact exactly replaces the temperature of the heat bath in the Einstein relation; 
it is different from the external temperature, but it can be shown to have all other 
properties defining a true temperature 14 . In particular, the measure of T^yn can 
be done using any version of the FDT relating the correlation and response of an 
observable ■ Subsequently, the violation of FDT has become a widely studied tool 
in the context of glassy dynamics ^| ED • 

The analogies between glassy thermal systems and non-thermal systems close 
to jamming ^S] have stimulated the investigations about the existence of dynamical 
or effective temperatures in non-thermal systems. In particular, the appearance of 
a dynamical temperature in models for dense, compacting granular media has been 
shown to arise from an Einstein-fike relation EOl EH E21 ESI ■ We wiU focus on this 
aspect in sectional 

Other works on dense granular matter have focused on another version of the 
Fluctuation-Dissipation Theorem which relates the energy fiuctuations to the heat 
capacity [241 1251 1^ . Moreover, another athermal system has been investigated by 
Ono et al. ■ With reference to a model of sheared foam, various possible definitions 
of effective temperature have been shown to coincide, in particular in the context of 
an Einstein-like relation and of energy fluctuations. 

On another side of the wide range of non-equilibrium systems, granular gases 
are very far from glassy systems. They are maintained in a dilute non-equilibrium 
steady-state in which the dissipation due to inelastic collisions between particles is 
compensated by a strong energy injection. However, the existence of an Einstein 
relation between diffusion and mobility of a tracer particle has also been investigated 
recently in this context EHIEniEniElj as will be developed in section O 

It should be remarked that up to now most of the work in this area has been 
theoretical and numerical and few experiments have been carried out in order to 
check the validity of the Fluctuation-Dissipation Theorem in real granular media. 
One example in this direction is discussed in section 01 concerning a vibro-fluidized 
medium in a liquid- like regime |32| . 
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3. Dense granular media 

3.1. A pedagogical 1 — d model 

In this section we consider a simple model which describes the evolution of a system 
of particles which hop on a lattice of k — 0, N stacked planes, as introduced in [21] 
and discussed in |34| . In particular the system represents an ensemble of particles 
which can move up or down in a system of N layers in such a way that their total 
number is conserved. We ignore the correlations among particles rearrangements and 
problems related to the mechanical stability of the system. The master equation for 
the density pk on a generic plane k, except for the k = plane, is given by: 

dtpk = (1 - Pk)D{pk)[pk-i -Pu + Pk+i ■ Pd)] + (1) 
- pfc[(l - Pk-i)D{pk-i)pd + (1 - Pk+i)D{pk+i)pu] 
where pd and pu (with pu +Pd = 1) represent the probabilities for the particles to move 
downwards or upwards, respectively, among the different planes. D(pk) represents a 
sort of mobility for the particles given by the probability that the particle could find 
enough space to move. Apart from other effects it mainly takes into account the 
geometrical effects of frustration, i.e. the fact that the packing prevents the free 
movement of the particles. With p„ and pd we can define the quantity x — Pu/Pd 
which quantifies the importance of gravity in the system. We can associate to x a 
temperature for the system given by T ~ 1 / log{l / x) . 

One interesting question to address is whether there exists a variational principle 
driving the relaxation phenomena in this system and in general in granular media. 
In other words one could ask whether some free-energy-like functional is minimized 
(Lyapunov functional) by the dynamical evolution. In this model, and more 

generally for evolution equations of the form 

dtPk = g{pk)[.f{pk-i)Pu+f{pk+i)Pd)]-.f{pk)[9{pk-i)Pd+9{pk+i)Pu]{'2) 

where / and g are generic functions for which we only require / > 0, 5 > 0, d/ /dp > 0, 
dg/dp < 0, it is possible to prove that such a functional, which decreases monotonically 
along the trajectories of the motion, indeed exists and is given by: 

oc 

F = Y,[j{x)kpk-S{pk)] , (3) 

fe=0 

with Sipk) = log g{p)/f{p)dp and -f{x) = log(l/a;). 

A deeper insight in the above mentioned phenomenology is obtained by 
considering the continuum limit for the model described by More precisely 

we consider a diffusive limit that consists in scaling the space variable as -, the 
time variable as and the drift term pd — Pu as e. Therefore x = ek, r — e'^k/2, 
Pd ^ Pu = and we consider the evolution of u{x) = p{k). We get the continuum 
limit by taking the Taylor expansion of the right hand side of jSJ around x = fee. For 
example p{k -I- 1) = u{x -I- e) = u{x) + edxU + ^e^dxxU + 0{e^). We get, formally, 

dru{x) = Pdxifg) + {gdxxf - fdxxg) + 0(e), (4) 
which, in the limit e ^ 0, gives 

dru{x) = I3dx{fg) + {gdxxf - fdxxg)- (5) 
This non-linear diffusion equation may be conveniently written in the following form 

/ dF\ 

dru = dx \ D{u)dx^ \ , (6) 
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where D = fg, ^ denotes the functional derivative of F with respect to u, 
F — /q°°(/3 u X — S{u))dx, and S' = log (^f^- Notice that the functional F decreases 
with the dynamics induced by One has, in fact, 

9rF ^ I dx^dru = / dx^d.. (D{u)d.^) (7) 



OU J OU \ OU J 

that, after an integration by parts, gives 



-Diu)\^-J <0. (8) 

Therefore there exists a "free energy" -like functional, F, for © which has exactly the 
same form of the functional defined for the discrete model (see Q). We can notice 
that, while the functional form of S and the value of /3 determine in a unique way 
the asymptotic state, they are not sufficient to determine the dynamical behavior of 
the system. In particular in order to know it one should know the functional form of 
D{p). 

What the analysis of this simple model suggests is the possibility of introducing, 
for non-thermal systems such as granular media, equilibrium concepts like free- 
energy, entropy and temperature. More precisely it is possible (in the case studied 
here) to predict the asymptotic state by means of the minimization of a suitable 
functional which can be constructed by entropic arguments. It is worth stressing how 
granular systems often exhibit memory and so the existence of a unique Lyapunov 
functional is not guaranteed in general. One could for example expect that several 
Lyapunov functionals are associated to different stationary states reached with 
different dynamical paths. 



3.2. Edwards' approach 

A very ambitious approach was put forward some years ago by S.F. Edwards and 
collaborators [321 1^ OHl OH] i by proposing for dense granular an equivalent of the 
microcanonical ensemble. The idea is to suggest that one could reproduce the 
observables attained dynamically by first measuring the density of the system, and 
then calculating the value of the observables in an ensemble consisting of all the 
'blocked' configurations at the measured density. The blocked configurations are 
defined as as those in which every grain is unable to move. 

This 'Edwards ensemble' leads naturally to the definition of an entropy SecIw, 
given by the logarithm of the number of blocked configurations of given volume, energy, 
etc., and its corresponding density SEdw = SecIw/N. Associated with this entropy 
are state variables such as 'compactivity' X^^^ — -^SEdwiV) and 'temperature' 

'^Edw = -^SEdw{E).^ 

That configurations with low mobility should be relevant in a jammed situation 
is rather evident, the strong hypothesis here is that the configurations reached 
dynamically are the typical ones of given energy and density. Had we restricted 
averages to blocked configurations having all macroscopic observables coinciding 
with the dynamical ones, the construction would exactly, and trivially, reproduce 
the dynamic results. The fact that conditioning averages to the observed energy 
and density suffices to give other dynamical observables (even if maybe only as an 
approximation) is highly non-trivial. 
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It turns out that the advances in glass theory mentioned in section [21 have in fact 
come to clarify and support such a hypothesis. Indeed, the dynamical temperature 
emerging from the Einstein-like relation between diffusion and mobility, despite its 
very different origin, matches exactly Edwards' ideas. One can indeed identify in 
mean-field models all the energy minima (the blocked configurations in a gradient 
descent dynamics), and calculate 1/TEdw as the derivative of the logarithm of their 
number with respect to the energy. An explicit computation shows that TecIw coincides 
with Tdyn obtained from the out of equilibrium dynamics of the models aging in contact 
with an almost zero temperature bath |4U[ 1411 I43L I45| . Moreover, given the 
energy E{t) at long times, the value of any other macroscopic observable is also given 
by the flat average over all blocked configurations of energy E{t). Within the same 
approximation, one can also treat systems that, like granular matter, present a non- 
linear friction and different kinds of energy input, and the conclusions remain the 
same 0^1 despite the fact that there is no thermal bath temperature. 

Edwards' scenario then happens to be correct within mean- field schemes and 
for very weak vibration or forcing. The problem that remains is to what extent it 
carries through to more realistic models. In the next subsections, we present the 
general methodology that has been used to explore this issue for some representative 
examples of models of granular compaction. 

3.3. General strategy to check Edwards' assumption: 

One possibility of making an assumption a la Edwards would be to consider a 
fast quench, and then propose that the configuration reached has the macroscopic 
properties of the typical blocked configurations. This would imply that the system 
stops at a density for which the number of blocked configurations is maximal. However, 
it turns out that generically the vast majority of the blocked configurations are much 
less compact than the one reached dynamically, even after abrupt quenches. 

One has therefore to give up trying to predict the dynamical energy or density by 
methods other than the dynamics itself. The strategy here is to quench the system to a 
situation of very weak tapping, shearing or thermal agitation. In this way, the system 
keeps compactifying, albeit at a very slow rate. Edwards' measure is then constructed 
as a flat measure over blocked configurations conditioned to have the energy and/or 
density of the dynamical situation. 

In order to check Edwards' hypothesis, we will show how Edwards' measure 
can be constructed in some representative (non mean-field) systems, together with 
the corresponding entropy and expectation values of observables. Moreover, we also 
construct what can be called Gibbs' measure by removing the constraint of sampling 
blocked configurations. Both measures are then compared with the observables 
obtained with an irreversible compaction dynamics. 

3.3.1. Models To illustrate the strategy to check Edwards' hypothesis we shall mainly 
refer to one specific finite-dimensional lattice model which has been shown to reproduce 
the complex phenomenology of the granular compaction. We focus on the so-called 
Kob- Andersen (KA) model 07|, which considers a three-dimensional lattice gas with 
at most one particle per site and otherwise purely dynamical constraints: a particle 
can move to a neighboring empty site, only if it has strictly less than m neighbors in 
the initial and in the final position. This model was first devised in the context of 
Mode-Coupling theories to reproduce the cage effect existing in supercooled liquids. 
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which produces at high density a very strong increase of the relaxation time|. Though 
very schematic, it has then been shown to reproduce rather weh several aspects of 
glasses 1321 ISO] and of granular compaction [ST] . The simplicity of its definition and 
the fact that it is non mean-field makes it a very good candidate to test Edwards' 
ideas: in fact, the triviality of its Gibbs measure allows to compare the numerical 
data obtained for the dynamics and for Edwards' measure with the analytic results 
for equilibrium. 

Completely similar results have been obtained in the framework of another class 
of non mean-field models, the so-called Tetris Model [6'6l 15] . Here the constraints are 
not purely dynamical, but related to the steric properties of the grains which undergo 
a geometrical (and hence dynamical) frustration. We refer the reader to HOI I^H HI] 
for details. 

3.3.2. Equilibrium measure Let us first consider the case of the Kob- Andersen model: 
the dynamic character of the rule guarantees that the equilibrium distribution is 
trivially simple since all the configurations of a given density are equally probable: 
the Hamiltonian is just since no static interaction exists. Density and chemical 
potential are related by p = 1/(1 + cxp(— /9/i)), and the exact equilibrium entropy 
density per particle reads 

Seg™i(p) = -plnp- (1 - p) ln(l - p) ^ ^^f^^ = (9) 

In this model, the temperature 1//3, which appears only as a factor of the chemical 
potential, is irrelevant, so that we can set it to one. Besides, the equilibrium structure 
factor is easily seen to be a constant, gequa{f) = P^- no correlations appear since the 
configurations are generated by putting particles at random on the lattice. It will 
therefore be easy, as already mentioned, to compare small deviations from gequa{T)i a 
notoriously difficult task to do in glassy systems. 

3.3.3. Edwards' measure Since Edwards' measure considers blocked configurations 
in which no particle is allowed to move, the crucial step to sample this measure is in 
fact to introduce an auxiliary model ^H]- the auxiliary energy Eaux is defined as the 
number of mobile particles^ where a particle is defined as mobile if it can be moved 
according to the dynamical rules of the original model. Edwards' measure is thus a flat 
sampling of the ground states {Eaux = 0) of this auxiliary model, which is obtained 
by an annealing procedure, at fixed density, of the auxiliary temperature Taux (and 
we write Paux = ^/Taux)- Note that the Monte-Carlo dynamics of the auxiliary model 
does not need to respect the constraints of the real model, so that efficient samplings 
with e.g. non-local moves can be obtained. 

From the measure of the auxiliary energy during the annealing, at given density 
p, EauxiPaux, p)i ouc cau computc the Edwards' entropy density defined by: 

SEdwip) = SauxiPaux = OO, p) = Sequilip)- / SauxiPaux , p)dPaux (10) 

where eaux{Paux,p) is the auxiliary Edwards' energy density and Sequuip) = 
SauxiPaux — 0, p). Fig. rcports the results for SEdw{p) as obtained from (|10|l 
compared with Sf.quii{p), for the Tetris (left panel) and the Kob- Andersen models 
(right panel), respectively. 

X For recent and interesting results about this model the reader is referred to 1481 
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Figure 1. Edwards' entropy density, s^dwip) and equilibrium entropy density, 
^Equil(p) for the Tetris (left) (see |20| for details) and the Kob-Andersen (right) 
models. 



The slope of the tangent to SEdw{p) for a generic p aUows to extract TEdwifi)- 
The natural definition for Edwards' temperature is 

_ 1 dsEdwjp) 

dp ' ^''^ 

which yields: 

„ , X dSequilip) ,dsEdw{p) 

J-EdwiP)- 1 / 1 • (.-L^j 

dp dp 

Similarly, the Edwards measure structure function, gEdwi^)^ is obtained as 

gsdwir) ^ lim gaux{r, Paux) ■ (13) 

3.3.4- Irreversible Compaction Dynamics The irreversible compaction dynamics is 
obtained by trying to increase the density of the system, starting from a low-density 
"equilibrium" situation. For the KA model, this can be done e.g. by increasing slowly 
the chemical potential on a given layer of a three-dimensional box: when the chemical 
potential becomes large enough, the dynamical constraint does not allow anymore the 
system to reach the desired density and slow compaction follows. 

During the compaction, we record the density p(t) and the density of mobile 
particles Pm{t). It is particularly interesting to notice that in the out-of-equilibrium 
configurations visited during the irreversible dynamics, the fraction of mobile 
particles is systematically smaller than the corresponding value in equilibrium. This 
suggests the possibility of distinguishing between equilibrium and out-of-equilibrium 
configurations by looking at the spatial organization of the particles in both cases. We 
have thus measured, during the compaction dynamics, the particle-particle correlation 
function at fixed density. 

The existence of an Einstein relation during the compaction dynamics is tested 
by the measure of the mobility of the particles 

a k=l •' 

obtained by the application of random forces to the particles between tw and -I- 1, 
and the mean square displacement 

1 ^ 

Bit^,t^ + = ^ E E + ^) - ^kitv.))^) , (15) 

a k=l 
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(N is the number of particles; a = 1, • • • , d runs over the spatial dimensions: d = 2 for 
Tetris, o? = 3 for KA; is measured in units of the bond size d of the square lattice). 
Indeed, the quantities x(tw,tw + t) and B{tyj,tw + t), at equilibrium, are linearly 
related (and actually depend only on t since time-translation invariance holds) by 

2x{t) = ^B{t), (16) 

where X is the so-called Fluctuation-Dissipation ratio (FDR) which is unitary in 
equilibrium. Any deviations from this linear law signals a violation of the Fluctuation- 
Dissipation Theorem (FDT). In particular, as mentioned in section 13 in many 
aging systems, and in particular in the KA model . violations from l|16() reduce 
to the occurrence of two regimes: a quasi-equilibrium regime with X — 1 (and 
time-translation invariance) for "short" time separations [t ^ t^), and the aging 
regime with a constant X < 1 for large time separations. This second slope is 
typically referred to as a dynamical temperature T^yn > T^q such that X — X^yn = 

3.3.5. Comparing different measures At this stage it is possible to compare 
equilibrium and Edwards' measures with the results of the out-of-equilibrium dynamics 
at large times. 

In Fig. 121 we plot the deviations of the dynamical particle-particle correlation 
functions from the uncorrelated value p^. In particular we compare {{gdyn{r) — P^)) 
obtained during the irreversible compaction with the corresponding functions obtained 
with the equilibrium and Edwards' measures. It is evident that the correlation 
function, as measured during the irreversible compaction dynamics, is significantly 
different from the one obtained with the equilibrium measure. On the other hand the 
correlation functions obtained with Edwards' measure are able to describe better what 
happens during the irreversible dynamics. In particular what is observed is that the 
correlation length seems to be smaller for configurations explored by the irreversible 
dynamics than in the equilibrium configurations. This aspect is captured by Edwards' 
measure which selects configurations with a reduced particle mobility. In practice one 
can summarize the problem as follows: given a certain density, one can arrange the 
particles in different ways. The different configurations obtained in this way differ in 
the particle mobility and this feature is reflected by the change in the particle-particle 
correlation properties. 

Another comparison can be performed with regard to the dynamical temperature 
Tdyn ISni- Figure 131 shows a plot of the mobility x(^i*tu) 'vs- the mean square 
displacement B(t,tu}) testing in the compaction data the existence of a dynamical 
temperature Tdyn [10]' The agreement between Tdyn and the Edwards temperature 
Tedun obtained from the blocked configurations as in Fig.Q] for the density at which 
the dynamical measurement were made, is clearly excellent. Further evidences in this 
direction have been obtained for the Tetris model |21| . 

To summarize, during the compaction, the system falls out of equilibrium at high 
density, and is therefore no more described by the equilibrium measure. It turns out 
that Edwards' measure, constructed by a flat sampling of the blocked configurations at 
the dynamically reached density, reproduces the physical quantities measured at large 
times, and in particular predicts the correct value for the dynamical temperature. 
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f r (distance) 



Figure 2. Left: KA model Dynamic structure function for the KA model, 
obtained in a very slow compression ("Dynamics"), and with the equilibrium 
("Gibbs") and Edwards' measure at the density reached dynamically. Right: 
Tetris model Comparison between the correlation functions obtained with the 
equilibrium measure, the Edwards measure {f3aux = 6) and the irreversible 
dynamics. In all cases the system is considered at a density of p ~ 0.58. 




B (t,t„) 

Figure 3. Einstein relation in the Kob-Andersen model: plot of the mobility 
Xitttm) VS. the mean-square displacement B{t,tm) (data shown as circles). The 
slope of the full straight line corresponds to the equilibrium temperature (T = 1), 
and the slope of the dashed one to Edwards' prescription obtained from figure [Tl 
at p{t^) = 0.848. 

3.4- Partial Conclusions 

It turns out that Edwards' measure, constructed by a flat sampling of the 
blocked configurations, is able to reproduce the physical quantities measured at 
large times. In particular, the connection of Edwards' temperature with the 
dynamical FDT temperature seems to be generally valid (though there are important 
counterexamples INDIES)- 

Apart from the first evidences reviewed here, which have lent credibility to 
Edwards' construction, various works have allowed a better comprehension of the 
validity of Edwards' approach and of its limitations. A comprehensive review of 
these approaches is beyond the scope of the present paper. An incomplete list would 
include results on 1-d models of particles or spins [21 UHl EH] : parking lot models |S5 j 
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3- d molecular dynamics simulations |23j . diluted spin-glasses and hard-spheres on 
lattice 1201, spin-glasses with tapping 

It is important to note that, in a case for which the explicit analytical computation 
of Edwards' measure and of the dynamical quantities is feasible, it has been shown 
that Edwards' construction of a flat measure on the blocked configurations is not 
exactly valid ■ In [2111 , the study of the low-temperature dynamics of the Random 
Field Ising model has also shown that the dynamically reached configurations have 
typically zero magnetization, in contrast to the configurations dominating Edwards' 
measure. While such studies show that Edwards' construction should probably be 
considered only as a first approximation, generalizations have been shown to yield 
better and better results for the prediction of observables obtained dynamically. These 
generalizations imply to use a flat average on blocked configurations, restricted by 
constraining more than just one variable to its value obtained dynamically |5HI I22j . 
The question arising concerns then the number of observables to constrain in order to 
obtain reliable predictions for the others. 

While such theoretical aspects remain interestingly open, a crucial question 
concerns the experimental validation of Edwards' ideas. In particular, the study of 
diffusion and mobility of different tracer particles within driven granular media would 
allow to confirm or disprove the theoretically predicted violation of FDT and the 
existence of dynamical temperatures §. Such experiments are actually in progress and 
a first set of results concerning the diffusion has been published in |58[ . 

4. Liquid-like granular media 

Context 

In this section we analyze whether a notion of temperature can be defined for a granular 
medium in a liquid-like regime, i.e. a regime where the medium is brought by vibration 
to a quasi-fluidized state. We shall mainly refer to some recent experimental results 
where a vibrate granular medium was sensed by means of a torsion pendulum. 

In the classical Brownian motion experiment, a "tagged" particle immersed in 
a liquid can be used as a thermometer to determine the temperature of the liquid 
itself. For this, one has to record the motion of the tagged particle, and data analysis, 
for example according to the Langevin formalism, gives the temperature. One may 
wonder whether a similar experiment, performed in a granular medium under suitable 
external vibrations (so that it looks very much like a liquid), could be used to determine 
a "granular-liquid temperature" . 

4- 2. Theoretical background 

Let us briefly review for clarity the behaviour of a torsion oscillator of moment of 
inertia / and elastic constant G immersed into an equilibrium liquid. Following the 
Langevin hypothesis |6U| , we suppose that the effect of this perturbing environment is 
split into two parts: a viscous friction force, proportional to the oscillator angular 
velocity, and a random, rapidly fluctuating force ^(t), which is an uncorrelated 
Gaussian white noise of zero mean and variance q. The oscillator angular position 

§ For other alternative experimental tests for the validity of Edwards' approach, inspired by the 
study of spin systems, see 1571 . 
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9 satisfies the Langevin equation: 

10{t) + a0{t) + G9{t) - ^{t) + C,^t{t) (17) 

where a is a friction coefficient and Cext denotes an external torque to which the 
system may also be submitted. When no external torque is applied {Cext = 0), a 
useful quantity that can be extracted from this equation is the power spectral density 
S, defined as twice the Fourier transform of the auto-correlation function {9{t)6{t')) 
((. . .) denotes the statistical average over the noise). Using the Wiener-Khintchine 
theorem for stationary processes, we get 

^('"^ ^ iH^S^^^ifTa^ ^^^^ 

where wq = \fGjl is the natural angular frequency of the oscillator. 

On the other hand, one can also focus on how the oscillator responds to an 
external torque Cext (t) ■ The quantity containing this information is the susceptibility 
x{t) (or linear response function), defined as 

9{t)^ fdt'x{t-t')Cext{t') (19) 



which implies that the external torque Cext{t) should be small enough for this linear 
approximation to be valid. From this definition and from the Langevin equation H17|) . 
we see that in the Fourier representation 

^^^^ " C^R " /(c^2_^2)+ic.^ (20) 

The real and imaginary parts of this complex function can be defined as x{^) = 
x'(w) — ix"{'-^), where in particular 

Comparing the power spectral density (|18|l and the imaginary part of the susceptibility 
H21|l. we now notice that these two very different concepts have similar expressions 
and are related by the simple relation 

X"(^) " a ^''> 
In a thermal system at equilibrium, using the equipartition of energy principle, 
the parameter q can be related to the bath temperature as g = 2akBT, thus giving 
the celebrated Fluctuation-Dissipation Theorem which states that 

S{lo)uj/x"H = 4fcBT. (23) 

Since a vibrated granular medium is not at equilibrium, there is no reason, in 
principle, to expect that a relation like (|23() should hold for such a system. We can 
expect that a granular medium can be found in quasi-stationary states but no ergodic 
principle can be invoked whatsoever. Nevertheless if a simple relation like H23|l were 
valid for a granular medium, at least in some regime, this would be a strong hint for 
the comprehension of the thermodynamic properties of such systems. 

The experiment we present was aimed precisely at the check of relation H23|l . 
The idea is to measure the noise power spectrum and the susceptibility while the 
granular medium is externally driven in the liquid-like state, and test whether there 
exists a Fluctuation-Dissipation-like relation. If so, this will give us a measure of a 
temperature-like parameter. 
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4-3. Experimental setup 

We use the following experimental setup |32| : a thin torsion oscillator is immersed at 
some depth into a granular medium made of millimeter-size glass beads, as shown in 
Fig. 0] (note the analogy with the situation described in for a system at thermal 
equilibrium). The beads are placed into a cylindrical container which is continuously 
vibrated vertically, with a high-frequency filtered white noise (cut off above 900 Hz 
and below 300 Hz in the experiments described) . We use this vibration mode to ensure 
a homogeneous agitation and avoid undesired effects such as pattern formation and 
convection rolls. Note that this type of white noise vibrations is not used in order 
to provide a random torque with white noise spectrum to the oscillator: actually, its 
motion is observed in a much lower frequency range (10-50 Hz) than the vibrations 
applied. 




Figure 4. Sketch of the torsion oscillator immersed into the "granular bath" . The 
granular medium, composed of glass beads of diameter 1.1 ±0.05 mm is placed in 
a cylindrical container of height 60 mm and diameter 94 mm. An accelerometer 
measures the intensity of the external perturbations, F. 



The vibration intensity is determined by an accelerometer fixed on the container, 
which measures the parameter F, defined as F^ — jA{f) df, where A{f) is the 
acceleration spectrum, normalized to the acceleration of gravity, and the integration 
is taken in the frequency range of about 1 Hz to 10 kHz. For sinusoidal vibrations, 
F = 1 is the "fluidization" threshold, above which a single grain starts to fly. Here, 
we typically use vibration intensities between F = 1 and 15. 

The oscillator angular position 9 is detected optically (see Fig. 01 . For 
susceptibility measurements, two external coils and a permanent magnet fixed on 
the oscillator axis allow to apply a sinusoidal torque Cext{t) = Ce sm{ujt). 



4-4- Results and discussion 

The analysis of angular deflection time-series 9{t) in the free mode {Cext — 0) provides 
the noise power spectral density S{llj), shown in Fig.Et. for different values of F. Then, 
with the oscillator in forced mode (with an externally applied torque Cext{t)), we 
measure the complex susceptibility x(w), whose modulus |x('-^)l is shown in Fig.Eb. 
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The amplitude of the external torque is small enough to be in the regime of linear 
response. 

Fitting these curves |x('^)| with the standard expression for the damped oscillator 
|x(w)| — [l"^ (wq — uj'^Y + a^oj'^]^^/'^ shows a good agreement, thus supporting the idea 
that the Langevin equation of motion (|17|l is a pertinent description of the oscillator 
linear response. This allows us to extract a granular friction coefficient a, or a granular 
viscosity /i oc a, found to be inversely proportional to the vibration intensity: a cx l/F 
(see inset of Fig. |S|d) . 




10 20 30 40 50 10 20 30 40 50 

f(Hz) /(Hz) 



Figure 5. Power spectral density (a) and modulus of the complex susceptibility 
(b) versus the frequency, / = a;/27r, for different vibration intensities F: from 
top to bottom, 11.6, 10.0, 7.3, 5.4, 3.7, 2.2, 1.5, 1.0. Results obtained with a 
conical probe of apex angle 120, covered with a single layer of glass beads. The 
immersion depth is about 11 mm. In (b), each curve is fitted (dashed line) with 
its corresponding equilibrium expression (see text for details). The moment of 
inertia of the oscillator is / = 1.5 X 10"^ kg-m-^, and the applied torque amplitude 
is Ce = 3.2 X 10~^ N-m. Inset: Friction coefficient a, obtained from the fit to the 
curves IxCi^)], versus F. The dotted line is a power law a oc l/F. 

We can also calculate the Fluctuation-Dissipation ratios S{uj)uj/{4:x"{uj)), which 
are shown in Fig. These ratios, even though not constant, are surprisingly 'flat', 
especially compared to what has been measured in other non-equilibrium thermal 
systems, such as in glycerol and laponite 

This reveals that the high-frequency driven agitation of the granular medium acts 
on the oscillator as a source of random torque with white spectrum, at least in the 
10-50 Hz range under consideration. Energy is thus injected at high frequency, and 
spreads into a low frequency white spectrum. 

Since those ratios do not exhibit a strong frequency dependence, this provides 
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support for the existence of a Fluctuation-Dissipation relation in off-equilibrium 
driven granular steady states. This relation can thus be used to define an effective 
temperature, Tf,ff. FigureEJ) shows the averaged Fluctuation-Dissipation ratio levels, 
that is, ksTeff, versus F (black curve). Fitting to a power law yields ksTgff cx F^ with 
p close to 2. This dependence suggests that, due to the complex dissipation processes 
between the grains, a fixed fraction of the energy input (vibrations) is effectively 
available as granular kinetic energy and is "sensed" by the oscillator. 




10 20 30 40 50 J 4 

/(Hz) r 



Figure 6. (a) Fluctuation-Dissipation ratios S(cj)tij/(4x"(i^)) versus the 
frequency / = uj/2tt for different vibration intensities (as in Fig. |FJ (b) 
Effective temperature versus vibration intensity F, as obtained from Fluctuation- 
Dissipation levels in (a) averaged between 10 and 50 Hz (black symbols), and from 
experiment using a conical probe with triple moment of inertia (red symbols). A 
power law fitted to the data gives T^ff oc F^, with p = 2.1. The dashed line has 
equation kgT^ff = 3.5 x lO^^^F^. 

In fact, we notice that the order of magnitude of the thermal energy /csTe//, 
as measured here, is consistent with realistic values of the mean kinetic energy 
per particle, as measured by grain-tracking methods |^ ES]- Thus, the effective 
temperature Tgff measured seems to be related to the granular temperature, as 
usually defined in granular gases. This is particularly interesting in view of the 
recent numerical experiments relating the granular temperature with a dynamical 
temperature obtained through FDT-like measures in granular gases (see section|Sland 

paEoi)- 

These results indicate that the use of macroscopic quantities, such as temperature- 
like or viscosity-like parameters, to describe the behavior of externally driven granular 
media is in first approximation possible. We can notice that, as shown in Fig.Eb, the 
temperature parameter obtained does not seem to depend on the moment of inertia of 
the probe used. On the other hand, we have also indications (not discussed in detail 
here) that several complications, related to the discrete and inhomogeneous nature 
of granular media and to the role of gravity, have to be taken into account. These 
mechanisms are triggered on as soon as the probe-granular medium interactions take 
place on length scales that are comparable to the grain dimensions, for instance when 
using non-smooth probe sections exhibiting features with length scales of the same 
order as the grain diameter, or close to the surface. . . The systematical study of these 
geometrical and grain-level effects will give us the possibility to establish empirical laws 
that may be used to predict granular behavior in practical situations, and possibly in 
simulations. 
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Let us note that the temperature measured here is not a priori the Edwards 
temperature. However, by decreasing the external vibrations we observe evidence 
of a glassy behavior ^21- In particular, the study of the power spectrum (hence of 
the diffusivity) reveals that for very weak external vibrations the granular medium 
exhibits a critical slowdown from the fluid state to a glassy state, where - as for 
supercooled liquids - the diffusivity approaches zero. In thermal systems undergoing 
a glass transition, the definition of a temperature is an open issue, and there are 
suggestions proposing that a "configurational temperature" , sharing similarities with 
the Edwards temperature, can be introduced. Thus, in the limit of very small external 
vibrations, it is possible that we are confronted with a similar problem in the granular 
media. 

5. Granular gases 

In this section the particular case of highly fluidized granular materials is taken into 
account. When a box full of grains is strongly shaken and the volume available is 
large enough, the assembly of grains behaves in a way similar to standard molecular 
gases. In the literature a large series of experiments has led to the definition of 
"granular gases", i.e. granular material in a gaseous state [SZlinHl- Because of the 
analogy with usual gases, the term "granular temperature" has naturally been defined 
as the average local kinetic energy per particle: this definition indeed coincides with 
the thermodynamic equilibrium temperature in the limit of elastic collision between 
particles. 

In experiments, it is possible to study kinetic as well as hydrodynamic observables 
and compare them with the results of statistical mechanics and hydrodynamics of 
molecular gases. The observed differences may be explained as a consequence of loss 
of energy in collisions between grains. Energy dissipation during a collision is due to 
irreversible transfer of energy from macroscopic energy to internal degrees of freedom 
and eventually to heat. Therefore it can be said that inelasticity is the main ingredient 
in the description of a granular material in the dilute regime, while frustration and 
excluded volume effects are negligible. We will show that, while inelastic collisions 
produce large deviations with respect to usual thermodynamics and hydrodynamics, 
linear response theory and Fluctuation-Dissipation relations are still valid |29[I3U[I5T] . 
provided that the "equilibrium temperature" is replaced by the granular temperature 
of each component of the gas. 

5.1. Fluctuation Dissipation relations for diluted granular systems 

We have studied Fluctuation-Dissipation relations using two sets of independent 
measurements, i.e. two choices of the pair response-correlation. 

Recipe I: The first one consists in the classical measure of mobility and diffusivity. 
The mean-square displacement 



behaves asymptotically as ~ 4D{t ~ t'). The mobility of a tracer particle can be 
measured by applying a constant and small drag force ^ — ^e^ to a given particle, 
labeled 0, for times t > t' (linearity of the response has been checked by changing the 




(24) 



What is the temperature of a granular medium? 



17 



amplitude of the perturbation). The perturbed particle will reach at large times a 
constant velocity fi, related to the response by 

X{t,t') = ^((ro(i) - ro(t')) • e^) « nt at large t. (25) 

At equilibrium Green-Kubo relations (Einstein relation in this case) predict /x = (3D, 
Tb = being the equilibrium temperature, so that x(^j*') = jB. 

Recipe II: Another totally independent way of checking FD relations is the 
following: once a steady-state has been reached, the system is perturbed impulsively at 
a given time to by a non-conservative force applied (non-uniformly) on every particle. 
The response is then monitored at later times. The force acting on particle i is 

F(r„f) =7«^(r„t) (26) 

with the properties Vx^7^0,V-^ = 0, where 7i is a particle dependent variable with 
randomly assigned ±1 values. A simple case is realized by a transverse perturbation 
^(r, f) = (0,£_cos{kxx)S(t)), where fcj. is compatible with the periodic boundary 
conditions, i.e. kx = 2Trnk/Lx with integer and the linear horizontal box size. 
The staggered response function R (i.e. the current induced at t by the perturbation 
at to), and the conjugated correlation C, 

' i 

C{t,to) = {J2Ut)Uto)coii{k4x,{t) - X,{to)]}) 

i 

are related, at equilibrium, by the FD relation R{t,to) = ^C{t,tQ), Tf, = 1/(3 being 
the equilibrium temperature. 

We have applied both recipes to a model of pure granular gas as well as to a 
model of granular binary mixture. We have also applied recipe I to a miich simplified 
model of pure granular gas (inelastic Maxwell model) where analytical calculations of 
diffusion and mobility can be straightforwardly obtained. 



5.2. Description of the models 

The simplest model of granular material in two dimensions is the hard disks gas with 
inelastic collisions. To counterbalance the loss of energy we enforce the stochastic 
forcing model, i.e. particles receive random acceleration as if they were in contact 
with a "heat bath" . Moreover, we may (or may not) provide a viscous drag to each 
particle. Viscosity has the role, in this model, of a regularizing force and allows for 
a better definition of elastic limit. It can be thought as the result of friction with 
external walls or with a fluid the gas is immersed in. However we will show that 
viscosity is not essential in this study and identical results can be obtained without it. 

We consider a volume V in dimension d = 2 containing = A^i + N2 inelastic 
hard disks, Ni and N2 being the number of particles in component 1 and 2 of the 
mixture, respectively. The disks have diameters a (identical for the two species) and 
masses mg. (where I < i < N and Si is the species index, 1 or 2, of particle i). In 
a collision between spheres i and j, characterized by the inelasticity parameter called 
coefficient of normal restitution as^sj , the pre-collisional velocity of particle i, Vj, is 
transformed into the post-coUisional velocity v- such that 

TTl 

<=^i- ^ IL (1 + ■ v,,)a (27) 
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where — Vj — Vj and a is the center to center unit vector from particle i to j 
{(^siSj = OisjSi so that the total linear momentum m^Vi + rrij'Vj is conserved). 

In between collisions, the particles are submitted to a random force in the form 
of an uncorrelated white noise (e.g. Gaussian) with the possible addition of a viscous 
term. The equation of motion for a particle is then 

m,^ -F, +R, -Cs.v, (28) 
at 

where is the force due to inelastic collisions, (st is the viscosity coefficient 
and {Ria{t)Rjp{t')) — ^I^5ij5ai35[t — t'), where Greek indexes refer to Cartesian 
coordinates. The granular temperature of species s is given by its mean kinetic energy 
Ts — ms{v'^) s / d where is an average restricted only to particles of species s. 

Pure system: when mi = m2 = m, Ci — C2 = Cj = = C and 
ail = Q!i2 = OL22 = Q;, the gas is monodisperse. This model has been extensively 
studied [SillTniinillllZniElESlEni. when C 7^ a "bath temperature" can be defined 
as Tf, = ^^/2C. This corresponds to the temperature of a gas obeying equation (|27|l 
with elastic collisions or without collisions. The same temperature can be observed if 
the viscosity is very high, i.e. when C 3> 1/tc where Tc is the mean free time between 
collisions. Here we recall that when a < 1 the gas still attains a stationary regime, 
but its granular temperature is smaller than Tf, and therefore the system is out of 
equilibrium. Moreover, the statistical properties of the gas are different from those of 
an elastic gas in contact with a thermal bath: mainly the velocity distribution is non- 
Gaussian with enhanced high-energy tails. The system is usually studied by means of 
molecular dynamics (with hard core interactions) or by means of numerical solutions 
of the associated Boltzmann equation. As we are interested in the dilute case, where 
Molecular Chaos is at work, we follow this second recipe, implementing the so-called 
Direct Simulation Monte Carlo 77 (DSMC) which is a numerical scheme that solves 
the Boltzmann equation for homogeneous or (spatially) non-homogeneous systems. 

Binary mixture: in the general binary mixture case, simulations as well as 
experiments and analytical calculations have shown that energy equipartition is 
broken, i.e. Ti 7^ T2. At the level of Boltzmann kinetic equation, the temperature 
ratio of a binary granular mixture subject to stochastic driving of the form given 
above has been obtained in jTS] for the case Cj. — and in 79 for C,si 7^ 0. In the 
case 7^ a bath temperature can still be defined as Ti, — ^^./2C^si- Note that in 
general and Cg. depend on mi and the correct elastic limit (i.e. equipartition) is 
recovered if and only if Tf, does not depend upon m^. In [7H| it has been shown that a 
model with oc ^rag^ and oc fairly reproduces experimental results for the 
temperature ratio T1/T2 measured in a gas of grains in a vertically vibrated box. It is 
also known that equipartition is not recovered even in the so-called tracer limit |8U| . 
i.e. in the case N2 = 1 and A'^i ^ 1. For binary mixtures we have implemented both 
molecular dynamics and DSMC, in order to study possible differences. 

Inelastic Maxwell Model: in analogy with elastic gases, some simplified model can 
be introduced. For instance, it has recently been proposed |H2 the inelastic analogue 
of a Maxwell gas [S^]. An elastic Maxwell gas is made by particles interacting through 
a special repulsive long range potential. The interest of such an interaction is that the 
corresponding collision frequency results strongly simplified: while for hard disks the 
collision frequency is proportional to the relative impact velocity 5, for particles of a 
Maxwell gas it becomes independent of 5. So the resulting Boltzmann equation takes 
a much simpler convolutive expression. The stochastic model of inelastic Maxwell 
gases is directly defined introducing a normal restitution coefficient smaller than one 
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in the kinetic equation of an elastic Maxwell gas. As in the elastic case, the interest 
for such operation is the achievement of interesting exact results. 

In the simplest case, i.e. in one dimension, the Boltzmann equation of an 
inelastic Maxwell gas reads dtP{v,t) = /3 J du P{u,t)P {Pv+{1 — P)u,t) —P{v,t), 
where /? = 2/(1 + a) and a is the restitution coefficient. At odds with the elastic case, 
its only stationary solution is the degenerate delta function, representing a system of 
particles at rest. Starting with a finite energy, the cooling of the gas has an easily 
computable rate, since the average squared velocity per particle decays exponentially 
voit)"^ = J v'^P{v,t)dv = voiOyexpi-Xt) with A = (1 - a^)/2. 

One of the interesting feature of this model, is that it admits a scaling solution of 
the form P{v,t) — i^j^fiv/vQ)^ although the scaling function / cannot be recovered 
through a simple analysis of the velocity momenta j81j . In fact, it can be shown that 
the scaling function that solves the equation, /(c) = f [i+esja i is a distribution whose 
moments higher than the second diverge. Note that, while the rate of dissipation 
depends on the restitution coefhcient, the scaling solution does not. This is peculiar 
of the one dimensional case. In higher dimension ,83, the algebraic tails of the solution 
are still present, but the exponent depends on the restitution coefficient (the tail is 
narrower and narrower reducing the inelasticity and the distribution becomes Gaussian 
in the elastic limit). The observation of scaling solutions for Maxwell gases is not new. 
However, for the elastic case, the solutions are not relevant, while in the inelastic case 
an initial distribution (e.g. an exponential or uniform) rapidly converges to the scaling 
solution. If we consider an initial distribution with a finite energy per particle and we 
define f{c,t) = vo{t)P{cvo,t), then the equation for the inelastic Maxwell model can 
be recasted into 

dtfic, t) + Xd, {cf{c, t))=pj dsfis, t)f{f3c + (1 - f3)s) - f{c) (29) 

This equation can be read as the Boltzmann equation of a new system, which is an 
inelastic Maxwell gas submitted to a special driving bath (it is often called "Gaussian 
thermostat") given by the term proportional to dc{cf{c,t)) [51]. Such a gas performs 
a stationary dynamics where the energy lost by inelastic collision is compensated by a 
negative viscosity term, which pushes the particles with a force proportional to their 
velocity. As we shall see in the following, it is possible to straightforwardly compute 
the mobility and the diffusion coefficient for a simple stochastic model governed by 
equation H29|) . This allows to explicitly check the validity of the Einstein relation for 
such a stationary non-equilibrium dynamics. 

5.3. Results 

Pure systems: in Fig. [TJ left frame, a parametric plot of mobility versus diffusion 
is displayed for several choices of parameters, showing the linear behavior analogous 
to Green-Kubo formulas (Einstein relation in this case). The same linear behavior 
is recovered plotting response versus perturbations in the case of an impulsive shear 
perturbation (recipe II experiment, see Fig. right frame). From the slope s of the 
observed linear behavior in the response-perturbation graph, one can get the effective 
temperature Te// = l/(4s). We always find Tejj = Tg, with Tg <Tb. 

Binary mixtures: in the case of binary mixtures, linearity or response- 
perturbation relations is again verified, in the mobility-diffusion experiment as well 
as in the current-shear perturbation experiment j3U| . Here we review just the former 
results, i.e. those for the Einstein relation. By successively using as test particle 
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Figure 7. Pure system. Left: parametric plot of x(*i*™) "^s. B{t,tw) for tlie 
numerical experiment witli recipe I (constant force applied on one particle) with 
r = 1, r = 0.8 and r = 0.7, with heating bath, and for different choices of the 
intensity S of the perturbation, using T;, = 0.1 ad r;, = 10, N = 500, Tc = 1, 
tm = 100. The results are obtained by averaging over 10000 realizations. Right: 
parametric plot of R{t — tw) vs. C{t — tw) for the numerical experiment with 
recipe II (impulsive shear perturbation) with r < 1, with heating bath, and for 
different choices of the wave number nj. of the perturbation. Tj, = 1 and t;, = 10, 
A'^ = 500, Tc = 1, S = 0.01, rifc = 8, with averages over 10000 realizations, using 
= 100. 



one particle of each species, one obtains the two responses xi and X2, and thus 
the mobihties fii and ^2- Two independent Einstein relations {^i = 2Di/Ti) are 
verified, by plotting Xi vs. Bi. In Fig. |Slwe show, as an example, the check of the 
validity of the Green-Kubo relations using DSMC in spatially homogeneous regime. 
All the experiments, performed varying the restitution coefficients and the masses of 
the two components, and with different models and algorithms (homogeneous and 
inhomogeneous, DSMC and MD) showed identical results, i.e the linearity of the 
response-perturbation relation with the effective temperature equal to the granular 
temperature of the perturbed species. 

In Fig. El an even more striking result is portrayed: the mobility-diffusion 
parametric graph is shown in the case of a single tracer with different properties 
with respect to a bulk gas (A^i — 500, A^2 = !)• In this case the tracer does 
not perturb significantly the bulk. However the temperature of the tracer is quite 
different from the bath temperature as well as from the gas temperature [80) . Again, 
the effective temperature of the tracer correspond to its temperature and not to the 
temperature of the bath or of the bulk. This is equivalent to say that a non-perturbing 
thermometer, used to measure temperature of a granular gas through Fluctuation- 
Response relations, would measure its own temperature and not the bulk temperature. 

Exact solution of the inelastic Maxwell model: 

Here we present a simple calculation on the inelastic one-dimensional Maxwell gas 
driven by a Gaussian bath where the (non-equilibrium) Einstein relation suggested by 
numerical simulation holds. 

Let us consider a gas of particles performing binary inelastic collision with a 
constant collision rate (i.e. independent of the precoUisional relative velocity) and 
submitted to a negative viscosity force, which drives the gas in a stationary state. 
The variation of the velocity of a particle after a time Ai is given by: 

, . , , , _ r Xv{t)At with prob. 1 — At , 

V[t+/\t)-v(t) - I _ l±a (^,(^) _ ^) ^ ^i^l^ p^Qb. At 
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Figure 8. Binary mixture, homogeneous DSMC: mobility vs. mean-square 
displacement; left: on = 0.3, ai2 = 0.5, 022 = 0.7, r?i2 = 3mi, Ti 0.2, 
T2 « 0.38; rieht: an = 0,9 = ac-c. = 0.9. mo. = 5mi . Ti ^ 0.035, T2 ^ 0.05. 
Sy: 2T2). 




Figure 9. Binary mixture, homogeeous DSMC: mobility Vs. diffusion of a single 
particle of mass mtracer in contact with N = 500 particles of mass m, immersed in 
a heat bath (i.e. random kicks plus viscosity). We use the following conventions: 
cttracer = Q:i2 and « = an. Only in case A) the tracer is also in contact 
with the external driving heat bath. A) mtracer = m, a = 0.9, cetracer = 0.4, 
Tg = 0.86, Tg"-'"=<='- = 0.70; B) mtracer = m, a = 0.9, atracer = 0.4, Tg = 0.86, 

j^tracer ^ o.46; C) mtracer = 7m, a = Otracer = 0.7, Tg = 0.74, T*'''"^'^'- = 0.60; 

D) mtracer = im, a = Otracer = 0.7, Tg = 0.74, Tg"'"'^'='' = 0.57. The solid line 
has slope J'*''"cer^ ^Yie dashed line has slope Tg. 



where u is the velocity of a generic colhding particle, distributed as f{u). The 
autocorrelation function A{ti,t2) = v{ti)v{t2) (with t2 > ti) can be computed 
using (|Sn|l . 

d 



-A(ti,t2) = lim v{ti) 

dt2 At^O 



V{t2 + At) ~ V{t2) 

At 



A 



l + a 



A{tl,t2), 



Since we are considering the stationary case, which is obtained with a negative viscous 
force exactly balancing the dissipation, i.e. A = (1 — a^)/4, this gives: 

A{ti,t2) ^Tg cxp - f ii^) (t2-ii) 



The diffusion coefficient can be computed via the autocorrelation function, obtaining: 



D = lim / dT'A{ti,ti + t') = — 

t^ccJq (1 



(31) 
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Now, to compute the mobility, we have to apply a small constant force to a particle, 

which modifies equation (|30|l so that ^^^^ = — '"^^^^ v{t) + F. The asymptotic velocity 
of the tracer yields 4F/(1 + a)^. This means that the mobility is 

4 _ D 

" + - 

and the (non-equilibrium) Einstein relation holds. 

5. ^. Summing up results about granular gases 

Approximated analytical results concerning Fluctuation-Dissipation relations in 
granular gases have been obtained for the case of mobility and diffusion of a tracer 
particle in homogeneous cooling granular gases (i.e. without external driving) in |86| 
and in driven granular gases in |31|. In the cooling case, two main ingredients spoil 
the usual Grecn-Kubo formula and lead to a strong reformulation of Fluctuation- 
Response relations: they arc the strong non-Gaussian behavior of velocity distribution 
(in Homogeneous Cooling State, HCS, the velocity p.d.f. has exponential tails), and 
the non-stationarity due to the thermal cooling of the gas. In the HCS, therefore, 
the Green-Kubo formula must be replaced by a non-linear formula which takes into 
account these two strong non-equilibrium effects. On the other side, the study of 
the tracer kinetic equation (Boltzmann-Lorentz equation) for the driven case leads 
to the conclusion that the only source of deviation from linearity in the Fluctuation- 
Response graph may be the velocity non-Gaussian behavior. However this is never 
very pronounced in granular gas driven by a homogeneous source (e.g. when grains are 
on a table which is vertically vibrated), so that deviations from linearity are negligible. 
The replacement of Tf, by Tg, comes quite naturally in the calculations. On the other 
side we have shown that in the inelastic Maxwell model, in Id, where the asymptotic 
velocity p.d.f. is known to have power law tails, the FD relations is recovered thanks 
to fortuitous balance of different terms. 

The general lesson learnt from simulations and analytical calculation, in the case 
of stationary dilute granular gas, is that FD relations are difficult to be violated when 
ergodicity is at work with mostly one characteristic time dominating the dynamics (i.e. 
the collision time). In such systems, the only cause of (always very small) deviations 
from usual FD relations is the non-Gaussian behavior of the velocity statistics, but 
even in particular cases where velocity pdf is strongly non-Gaussian Green-Kubo 
formulas may possibly work. The fact that the granular temperature of the measured 
tracer is the effective temperature is a quite obvious result if the original derivation 
of Green-Kubo relations is followed, as there the effective temperature appears to be 
simply (w^). 

6. Conclusions 

In this paper we have briefly reviewed the different approaches which have been 
followed in the last few years in order to define a notion of temperature for granular 
media. This question is a non-trivial one for such systems where the usual notion 
of termodynamic temperature, related to thermal agitation, does not play any 
obvious role. In these systems the very possibility of consistently constructing a 
thermodynamics is doubtful due to the fact that energy is lost through internal friction, 
and gained by non-thermal sources such as tapping or shearing. The dynamical 
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equations, whenever one could be able to write them down explicitily, do not leave the 
microcanonical or any other known ensemble invariant. Moreover these systems could 
never be considered at equilibrium and even the existence of stationary states is not 
always guaranteed. For instance often for dense granular media, just as in the case of 
aging glasses, a stationary state cannot be reached on experimental time scales. 

Despite all these difhculties in the last few years there have been several 
contributions which, though not yet completely satisfactory, are interesting because 
they have opened a new perspective which is worth pursuing in the future. 

One of the key concept has been the notion of dynamic (or effective) temperature, 
as defined in the framework of the Fluctuation-Dissipation Theorem. Following the 
remarkable work done on glassy systems where the partial violation of the Fluctuation- 
Dissipation Theorem has been put in relation with the existence of a so-called 
dynamical temperature (describing the slow structural rearrangements of the system), 
a lot of work has been done along the same lines in the framework of granular media. 

Another key contribution is due to S.F. Edwards who put forward a very 
ambitious approach to define a granular "ensemble" by looking at the so-called blocked 
(or jammed) configurations. 

In this paper we have tried to sum up (in a very partial and maybe subjective 
way) all these efforts, classifying them with respect to the different regimes a granular 
medium can be found in: the glassy regime, the liquid-like and the granular gas one. 

For the glassy-like regime in particular, we have reviewed Edwards' approach 
and described a possible path to check its validity for two non-mean field models: 
the Kob-Andersen model and the Tetris model. From this and other studies it 
turns out that the notion of Edwards compactivity seems to be closely related to 
that of dynamical temperature. A somewhat surprising but very interesting result, 
especially because it opens the way to experimental checks of Edwards' hypothesis. 
For the liquid-like regime we have reported about recent experimental results where 
an unusual "thermometer" , in the specific case a torsion pendulum, has been used to 
test the soimdness of the temperature concept in a continuously shaken container of 
tiny beads. Also in this case the temperature has been defined in the framework of 
the Fluctuation-Dissipation Theorem and its value seems to be consistent with values 
of the so-called granular temperature, defined in terms of the velocity fluctuations. 
The relation of this temperature with the one defined in the glassy regime is an open 
problem. Finally for the granular gas regime, we have reported about the validity 
of the Fluctuation-Dissipation Theorem. It turns out that, while inelastic collisions 
produce large deviations with respect to usual thermodynamics and hydrodynamics, 
linear response theory and Fluctuation-Dissipation relations are still valid provided 
that the "equilibrium temperature" is replaced by the granular temperature of each 
component of the gas. There exist deviations with respect to the usual Green-Kubo 
formula which are due to non-Gaussian behavior of velocity distribution (in particular 
for cooling granular gases) and the non-stationarity due to the thermal cooling of the 
gas. These deviations are indeed very small for driven granular gases. 

The picture emerging is still partial, even in each specific regime. It would be 
important to reinforce in the next years the experimental research in order to check the 
theoretical predictions and try to bridge some links between the different regimes, even 
though we expect that the level of universality, for these non-equilibrium systems, is 
very low. The two extreme regimes seem the most lacking. In particular in the 
glassy regime it would be important to have some experimental check of Edwards's 
hypothesis. On the other hand for granular gases there are already many predictions 
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which only call for an experimental check. 
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